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Recent progress in the theoretical description of inclusive B — > X u l~v decays in the shape-function region is 
reported. Finite moments of the shape function are related to HQET parameters. Event fractions for several 
experimental cuts are presented, with a particular emphasis on the hadronic variable P+ — Eh — \Ph\- The 
aim of this talk is to introduce the P+ spectrum, to compare it to the hadronic invariant mass spectrum and 
the charged-lepton energy spectrum, and to study the prospect of evaluating \V u b\ in the presence of a large 
background from B — > X c decays. 



1. INTRODUCTION 

Studies of both exclusive and inclusive semilep- 
tonic decays of the B meson can be used to ex- 
tract the magnitude of the CKM matrix element 
\Vub\- Both methods require some input from 
theory. Because the exclusive modes suffer from 
large form-factor uncertainties, the determination 
from inclusive B — > X u l~v decays are theoreti- 
cally favored. The overall relative uncertainty on 
\V u b\ measurements is currently about 15% [1], 
and there is hope for a significant reduction of 
the dominant theoretical errors in the future. The 
main problem in the computation of the relevant 
quantities is that it requires a framework which 
includes a systematic treatment of both pertur- 
bative corrections (including Sudakov resumma- 
tion) and non-perturbative effects (shape func- 
tion). Much progress in QCD-Factorization of 
inclusive B decays has been made in the past sev- 
eral months [2-4], and the aforementioned prob- 
lem is now understood using a sophisticated effec- 
tive field theory machinery. The main ingredients 
are Soft-Collinear Effective Theory (SCET) [5-7], 
Heavy-Quark Effective Theory (HQET) [8], and 
Renormalization-Group (RG) evolution of bilocal 
operators on the light cone [9] . 

Let us begin by discussing the kinematic setup 
used in the description of a typical B — ► X u \~v 
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Figure 1. Hadronic phase space in the light-cone 
variables P+ and P_ . The dark gray area is free 
of B — > X c background events, which occur in 
the light gray wedge above. For an explanation 
of the dashed and dotted curve, see the text. 



event. The momentum of the hadronic X u system 
is denoted by P H , and we work in the rest frame of 
the B meson. It is useful to describe the hadronic 
phase space in the variables P + = Eh — \Ph\ and 
P_ = Eh + \Ph\, in which it takes the simple 
form of a triangle (neglecting the pion mass) 

< P+ < P- < M B ■ (1) 
The phase space is depicted in Figure 1. The use 
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of the light-cone variables P + and P_ makes it 
easy to interpret an event (a point) in this pic- 
ture: the region of low recoil is along the diag- 
onal, while an event with a maximally recoiling 
X u is located on the right. Events with final-state 
hadronic invariant masses larger than Mo are lo- 
cated in the light-gray wedge in the upper part 
of the Figure, since M x = P+P-. The dark-gray 
region underneath is free of the B — > X c back- 
ground, which is about 60 times larger than the 
B — > X u signal. 

There are many possibilities to eliminate the 
background by cutting on kinematic variables. 
Let us briefly mention a few of them. Clearly 
a cut on M x < M D is the "ideal" separator, de- 
noted as a solid line in Figure 1 . The dashed hor- 
izontal line is located at P+ = M d /Mb, which 
marks an alternative way of eliminating charm 
events. A cut on the charged-lepton energy Ei > 
(M§ - M D )/2M B samples the same hadronic 
phase space as a cut on P + < M d /Mb- However, 
it contains far fewer events. Finally, a cut on the 
dilepton momentum squares q 2 < (Mb — Mo) 2 
leads to the area underneath the dotted line. 

The phase space is most densely populated 
in the "shape-function region" of large P_ ~ 
0{M B ) and small P+ ~ 0(A QC d)- In order to 
determine \V u b\ from this inclusive decay mode it 
is necessary to know what fractions of events sur- 
vive the various experimental cuts. In this talk we 
report on progress to answer this question. In ad- 
dition, we advertise cutting on P + as an efficient 
method for future \V u b\ determinations. To this 
end we study the advantages and disadvantages of 
this method over comparable ones, in particular 
the "ideal separator" Mx- 

2. FACTORIZATION OF THE DIFFER- 
ENTIAL DECAY RATE 

Recently, a systematic framework has been de- 
veloped [2-4] that enables us to compute the dif- 
ferential decay rates in the shape-function region. 
A series of matching calculations QCD — > SCET 
— > HQET is necessary to disentangle physics at 
the three different energy scales mj, -J mb Aqcd, 
and Aqcd ■ The general methodology is schemat- 
ically visualized in Figure 2. 
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Figure 2. Contributions from different energy 
scales are factorized by matching onto Effective 
Field Theories: QCD -> SCET -> HQET. 



The idea is as follows: Under the assump- 
tion of quark-hadron duality, the inclusive de- 
cay rate can be calculated in QCD using the 
optical theorem. The intermediate line denotes 
the energetic u-quark propagator. At a hard 
scale fa ~ m bi the amplitude is matched onto 
SCET, which correctly describes the infra-red de- 
grees of freedom below this scale. Physics ef- 
fects from scales above fa are contained in a 
hard coefficient function H(fa, m-b)- Since the it- 
quark momentum is parametrically off-shell by an 
amount of 0(m& Aqcd), the corresponding prop- 
agator can be integrated out at an intermediate 
matching scale fa ~ \J ti&Aqcd, leading to a per- 
turbatively calculable jet function J and a low- 
energy description in HQET. The matrix element 
of the resulting leading power operator defines the 
hadronic shape function, which cannot be com- 
puted using analytic techniques. Large (Sudakov) 
logarithms are resummed when evolving the hard 
function from the hard scale /j,h down to the in- 
termediate scale fa. Finally the decay amplitude 
is expressed as a convolution integral of perturba- 
tively calculable functions and the shape function 
renormalized at fa. 



3 



Event fractions and spectra are derived by per- 
forming the necessary phase-space integrations, 
before integrating over the shape function. In 
this way, the results are model independent and 
once again given as convolution integrals over the 
shape function. 

3. PROPERTIES OF THE SHAPE 
FUNCTION 

The shape function is a non-perturbative struc- 
ture function that encodes the Fermi motion of 
the heavy quark inside the B meson [10]. Al- 
though the functional form of it cannot be derived 
using analytical techniques, it is nevertheless pos- 
sible to study other properties of it using pertur- 
bation theory. In particular, the dependence on 
the renormalization scale /z can be reliably com- 
puted as long as [i is much larger than Aqcd • At 
leading order in renormalization-group improved 
perturbation theory it is given by the simple for- 
mula (for jUj > /Jo) 



S(u>, Hi 



S(u>',H)) 



(2) 



where r\ = (16/25) ha.(a s (no)/at s {Hi)), and the ex- 
ponent Vs vanishes in the limit fii — > fj,o- (The 
hatted notation denotes that the shape function 
is defined in a scheme-independent way and has 
support for cD £ [0,oo[. For details, see Ref. [3].) 
As a consequence of (2), radiative corrections 
build up a tail that vanishes slower than u~ l for 
w — > oo. This means that moments of the shape 
function (including its norm) are UV-divergent. 
This is, however, not an obstacle in practice, 
where the shape function is only needed over a 
finite interval. It is then natural to define mo- 
ments accordingly over a finite integration do- 
main, by means of a cutoff. Expanding them 
in a local OPE allows us to extract information 
about the shape function from the measurement 
of HQET parameters such as A and fi^. The 
shape-function mass scheme [3] is convenient for 
such calculations and free of renormalon ambigu- 
ities. We stress, however, that information of A 
and /ij in any low-scale subtracted mass scheme 
is useful, since there exist perturbatively calcula- 
ble relations between these mass definitions. 



A somewhat surprising observation is made by 
inspecting the behaviour of the moments under 
variation of the cutoff. It turns out that their 
values decrease for larger cutoffs. It follows that 
the tail of the shape function must be negative! 
In fact, it is possible to derive the explicit form of 
the shape function itself for large values of w > 
Aqcd- We find that the asymptotic tail is given 
in the MS scheme as [3] 



<i>A QCD C_pa s (/z) 
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(3) 



and one has to drop the common interpretation of 
the renormalized shape function as a probability 
distribution. 

3.1. MODEL SHAPE FUNCTIONS 

In the region uj ~ Aqcd, where it is rele- 
vant for phcnomenological applications, the shape 
function cannot be computed within perturbation 
theory In the following we will use a simple 
model that is consistent with all analytic con- 
straints mentioned above. The model is con- 
structed in such a way that its first few moments 
are consistent with the experimental values A = 
(0.63 ± 0.07) GeV and /4 = (0.27 ± 0.07) GeV 2 
at a reference scale of /ij = 1.5 GeV in the 
shape-function scheme (see Ref. [3] and references 
therein). Below we will use nine different func- 
tions (corresponding to the nine different pairs of 
values (A, ^ ) when varied within their errors) to 
account for our ignorance of the functional form 




u [GeV] 



& [GeV] 



Figure 3. Model shape function under variation 
of the numerical values for A and ^ within their 
errors. Correlated variations are shown on the 
left, anti-correlated on the right. 
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of the shape function. Figure 3 is intended to 
give the reader an idea of such variations. We 
stress, however, that a hypothetical knowledge of 
the first few moments to arbitrary precision would 
not determine the functional form of the shape 
function. 

4. EVENT FRACTIONS WITH CUTS 

ON P + OR M x 

With the methodology outlined in Section 2 
we find the following expression for the frac- 
tion of events that survive a cut P + < Ap, 
i. c. that are located below a horizontal line in 
the phase-space picture in Figure 1. At leading 
power in Aqcd/Mb and at next-to-leading order 
in rcnormalization-group improved perturbation 
theory, it is given as a weighted integral over the 
shape function: 

Fp(Ap) = T{a) e VH ^ h '^ x 

J du)S(u, Hi)Wp(u, A P ,a) . (4) 
o 

Large Sudakov logarithms are resummed in the 
expression T(a) e Vll ^ h '^ (where a = 16/25 x 
\n(a s (fj,i)/a s (fj,h)) ~ 0.3 for typical choices of 
Hh,^), an d the weight function Wp is perturba- 
tively calculable [3] and starts as Wp = 1+0 (a s ). 

When cutting on the hadronic invariant mass 
M x < s instead, the corresponding event frac- 
tion is given by 

F M {s ) = F P {s /M B )+T(a)e VH ^^ x 

J duS(u, Mi)Wtriangle(w) . (5) 

The two terms in the above expression can be 
easily interpreted using the phase-space picture 
in Figure 1. For an optimal cut so = M|, the 
allowed region (dark gray) splits into two dis- 
tinct areas, separated by the dashed line. Events 
below that line are included in the count by 
the first term in (5). The second term counts 
events in the triangle-shaped region above the 
dashed line. To account for the tip of this tri- 
angle (P + = P- = y / so) is problematic in the 



aforementioned framework, because the assump- 
tion P + <C P- and the resulting power-expansion 
rules are no longer valid. However, only few 
events are located near the tip, which is reflected 
in the finding that the corresponding contribution 
to F M scales like (Aqcd/Mb)^ 01 )/ 2 (assuming 
so ~ Aqcd-Mb)- Therefore one obtains a well de- 
fined leading power expression by setting the up- 
per limit of integration from ^so to oo [3]. It is 
unclear, however, how to treat power corrections 
to Fm systematically in the future improvement 
of this result. 

4.1. MODEL-INDEPENDENT 
RELATIONS 

QCD-Factorization is a powerful tool because 
the non-perturbative structure functions entering 
the theoretical description are universal, process- 
independent quantities. The shape function 
needed in this discussion on B — > X u l~v de- 
cays is also the only leading-power structure func- 
tion entering the calculation of the important 
B — > X s "f decay rate. In fact, at leading power 
and leading order in a s , the shape function is 
identical to the photon spectrum. This implies 
that it is possible to construct relations between 
the photon spectrum and the event distribution 
functions in B — > X u l~v decays in which the 
shape function has been eliminated. An example 
of such a relation is [4] 

Mb. 

2 

F P (A)= J dE^^w^Ry). (6) 

M B~* 1+ Qs ln... 

The weight function w is perturbativcly calcula- 
ble once our formalism has been applied to calcu- 
late the B —* X s 7 photon spectrum. Note that 
a similar relation to the hadronic invariant mass 
distribution Fm would be possible, but more com- 
plicated. Because of the contribution from the 
triangle region, the photon spectrum would be 
needed over a larger window and beyond the re- 
gion where it is experimentally known with ac- 
ceptable precision. 

4.2. MODEL PREDICTIONS 

We use the nine model shape functions dis- 
cussed in Section 3.1 to predict the event fractions 
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A P [GeV] A s [GeV] 

Figure 4. Model Predictions for F P (A P ) and F M (A S ) with A s = s /M B . Each gray band corresponds 
to a fixed value of A(/ij,/ij), while the value of is varied within its error bars. The solid lines 

on the right in each graph denote the inclusive B — > X c background, while the dashed lines indicate 
exclusive B — > Z)W decays. 



Fp and Fm and demonstrate their sensitivity to 
the shape function. In principle, one also has to 
vary the functional form of the shape function 
without violating the constraints for its first few 
moments. In the end we find, however, that such 
variations are already covered in the sample of 
the nine models [3]. 

In Figure 4 our results are depicted as gray 
bands. Each of the bands represents the results 
using shape function models that correspond to 
a fixed value of A, while the numerical value of 
/ij is varied within its error bars. We draw three 
bands for three values of A: 0.56 GeV, 0.63 GeV, 
and 0.70 GeV, in accordance with our discussion 
in Section 3.1. The calculations are valid for Ap 
and A s = sq/Mb of order Aqcd- If these val- 
ues are too small, i. e. parametrically of order 
Aq CD /Mp, w 50 MeV, the assumption of quark- 
hadron duality is no longer justified. If, on the 
other hand, Ap ;S are larger than about 0.8 GeV, 
the collinear expansion used in the calculation 
breaks down. 

Our findings arc summarized in Table 1 be- 
low. For "ideal" cuts, i. c. where the charm 
background starts, A P , S = M^/Mb ss 660 MeV. 
The results show that cutting on the hadronic in- 
variant mass Mx < Mp> is rather insensitive to 
shape-function effects. However, due to detector 



resolution effects, it is typically necessary to lower 
the cut and move away from the point where the 
charm background starts. In that case, the uncer- 
tainties due to our ignorance of the shape function 
become quickly larger, and are of comparable size 
to the uncertainties in an "ideal" P + cut. It seems 
a clear advantage for cutting on hadronic invari- 
ant mass. However, we will argue below that it 
might not be necessary to move too far away from 
the charm background when cutting on P + . In 
that case the hadronic invariant mass cut does 
not offer that advantage anymore and both meth- 
ods give high efficiencies of about 70-80% with a 
relative uncertainty of roughly 10%. 



Table 1 

Shape-function uncertainties for typical cuts on 
the hadronic variables Mx and P + . 



Cut 


Efficiency 


M x < M D 
M x < (1.7 GeV) 
M x < (1.55 GeV) 


(81.4+^)% 
(78.2+^)% 
(72.7±g;|)% 


P+ < H| = 0.66 GeV 


(79.6±|"1)% 


P+ < 0.55 GeV 


(69.01?^)% 
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5. CUTTING ON THE CHARGED LEP- 
TON ENERGY 

A prominent alternative way to discriminate 
the charm background is to cut on the charged- 
lepton energy. This is experimentally favored be- 
cause it does not require the reconstruction of the 
neutrino momentum. (Both cutting on P+ or Mx 
need a neutrino reconstruction.) Unfortunately, 
the method is not favored by theorists, for several 
reasons. The most obvious one is that far fewer 
events survive such a cut. In our prediction this 
is reflected by the fact that the weight function is 
suppressed by a power of AqcdACb : 

Fe(Ae) = T(a) e VH ^ h '^ x 

f dw S(w, m) 2( .^~^ [1 + 0(a s )] , (7) 
J M B - lu 

o 

where = Mb—2Eq and Eq is the lower limit of 
allowed lepton energy. Because of the rather low 
efficiency, this method is more prone to uncertain- 
ties from other effects, e.g. weak annihilation [11]. 

Numerical predictions using shape-function 
models are visualized in Figure 5 and summa- 
rized in Table 2. For orientation, the beginning 
of the charm background at Ab = 660 MeV 
is marked with an arrow in the Figure. Evi- 
dently, the relative shape-function uncertainties 
are large, and do not show a "focusing effect" 
as in the case of P+ or Mx cuts. Because of 
the large shape-function sensitivity and the over- 
all small efficiency, the extraction of \V u b\ from 



Table 2 

Predictions and shape-function uncertainties for 
a cut on the charged lepton energy. For com- 
parison, the results using the DeFazio/Neubert 
(DFN) model [12] are also given. 



cut E 


A E 


DFN 


BLNP 


[GeV] 


[GeV] 


[%] 


[%] 


Ei > 2.31 
Ei > 2.2 
Ei > 2.1 


0.66 
0.88 
1.08 


7.9+lj 
14.4+ 4 / 4 
20.9+i? 


99 9+3.2 

3 

31.7+!;$ 




0.2 0.4 0.6 0.8 1 



A E = M B - 2E [GeV] 

Figure 5. Model Predictions for the event frac- 
tion Fe- The bands are explained in Figure 4, 
and the arrow indicates the start of the charm 
background. 



the charged-lepton energy endpoint region is the- 
oretically disfavored. 

An interesting observation can be made when 
comparing the third column in Table 2 with our 
results (BLNP = Bosch, Lange, Neubert, Paz [3]) 
in the fourth column. The former are predic- 
tions using the familiar DeFazio/Neubert model 
(DFN) of Ref. [12], where a simple replacement 
rule is used to implement the shape function into 
the results of parton- level calculations. The cru- 
cial difference is that our framework allows for 
a clean separation of short- and long-distance 
physics effects. In order to demonstrate this, con- 
sider the one-loop order result of the triple differ- 
ential decay rate [2,3,12]. The leading logarithm 
comes with a coefficient of —4. In the frame- 
work of factorization, this coefficient will split 
into —4 = —8 + 4. While the leading logarithm of 
the jet function comes with the coefficient +4, the 
—8 will be absorbed in the "partonic" shape func- 
tion. The simple replacement rule of the DFN 
model does not capture this effect correctly. In 
addition to this important observation, the fac- 
torization framework allows for a systematic Su- 
dakov resummation, which was not performed in 
the work of [12]. 

The numerical impact of these improvements 
over our previous understanding is that the pre- 
dictions for the event fractions in the charged- 
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lepton energy endpoint have increased quite no- 
ticeably. This implies in turn that experimen- 
tal \V u b\ determinations from the charged- lepton 
energy endpoint region might need to be cor- 
rected to lower values. In general, our analysis for 
B — > X u l~v decay spectra in the shape-function 
region suggests that more events than previously 
anticipated are located in the shape-function re- 
gion. It will be most interesting to study how 
power corrections will affect this conclusion. 

6. CHARM BACKGROUND 

The P+ spectrum has a considerable advan- 
tage over the hadronic-mass spectrum when con- 
sidering the charm background. It is straight- 
forward to visualize the area in phase space that 
is populated by inclusive B — > X c decays. The 
OPE prediction uses quark-hadron duality and 
the fact that p+p~ > m 2 in the parton picture. 
Here, p± = P± — A are the light-cone compo- 
nents of the parton momentum, and m c is the 
charm quark mass. The dotted line in Figure 6 
marks the threshold p+p- = p 2 = m 2 c . At tree- 
level, all events are located on that line. Radia- 
tive corrections smear the event distribution into 
the area above it. In this sense, the tree-level 
scenario serves as an upper bound of the inclu- 
sive charm background. Note that the dotted line 
touches the tip of the triangle, P + = P_ = Md 
and extends to the right while always staying 
above the exclusive B — > D solid line. It fol- 
lows that the variable P+ is bound to be always 
greater than approximately 960 MeV in the par- 
ton model. Therefore an ideal cut on P + = 660 
MeV has the advantage of a comfortable "buffer 
zone" to the inclusive background. 

The gap between the charm-free area (dark 
gray) and the inclusive charm background (above 
the dotted line) is filled with exclusive final states 
such as B -> D, B -> D* , B -> Dtt, etc. 

In Figure 4 the results for a tree-level calcula- 
tion [4] of the inclusive and exclusive event frac- 
tions are drawn as solid and short-dashed lines, 
respectively. Note that the event fractions are 
normalized to unity; however, the reader should 
keep in mind that the B — > X c background is 
about 60 times larger than the B — > X u signal. 
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Figure 6. Phase-space picture for the charm back- 
ground. The exclusive B — > D events are lo- 
cated on the lower solid line. The OPE prediction 
for inclusive B — ► X c events are located above 
the dashed line. The gap is filled with exclusive 
events B — > D* , B — > Dir, etc. 



Consider now what happens when the quanti- 
ties Ap jS approach the charm threshold near 660 
MeV. The charmed hadronic invariant mass spec- 
trum receives step increments when the D and 
D* phase space becomes available. At the same 
time, the inclusive rates start near 660 MeV and 
reaches unity at 960 MeV in the tree approxima- 
tion. Detector resolution effects smear the charm 
events into the region below 660 MeV, and force 
us to move away from the ideal cut, which leads 
to enhanced shape-function uncertainties in the 
B — > X u event fractions (see Table 1). 

The scenario for the P + spectrum is different. 
The short-dashed line in the left part of Figure 4 
also starts near 660 MeV, but has a smooth onset. 
The same applies when the D* phase space be- 
comes available, leading to a (tiny) kink in the 
line. However, up until this point, the charm 
background consists exclusively of B — > D events, 
since the inclusive events start much later at 960 
MeV. The knowledge of the precise nature of the 
background ought to help experimentally to stay 
closer to the ideal cut, because it can be mod- 
eled with greater precision. We believe that this 
fact, together with the cleanliness of the theoreti- 
cal description, make the P + method an excellent 
high-efficiency candidate for a determination of 
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| V u f, |, and we urge the experimental community 
to perform a measurement of the P + spectrum. 

7. CONCLUSION 

In this talk we have presented recent advances 
in the understanding of inclusive B decays in the 
shape-function region. Using effective field the- 
ory techniques, a systematic framework was de- 
veloped, where a simple OPE cannot be used 
due to soft shape function effects. New results 
for event distributions of several important kine- 
matic quantities, including the hadronic invariant 
mass and the charged-lepton energy, have been 
presented. We advertised a new method for the 
extraction of \V u b\ based on the P + spectrum. 
We stress, however, that this should not be un- 
derstood as a discouragement of pursuing other 
methods, but rather as an encouragement to con- 
sider the P + method as a powerful alternative. 

The calculations have been performed in the 
heavy-quark limit and at next-to-leading order in 
RG-improved perturbation theory. Our predic- 
tion for the fraction of events with the optimal 
cut P + < Ml/ 'M B is [4] 

F P = (79.6 ±10.8 ±6.2 ±8.0)%, (8) 

where the errors represent the sensitivity to the 
shape function, an estimate of 0{a 2 s ) contribu- 
tions, and power corrections, respectively. Each 
of these uncertainties can be reduced in the future 
within the general framework presented in [2-4] 
and outlined in Section 2. Much can be learned 
about the shape function from the B — > X s j 
photon spectrum, once the formalism has been 
applied. The calculation of 0{a 2 s ) corrections 
would be a major effort, but is feasible in the 
future. The important question of power correc- 
tions are currently investigated by several groups 
in the theory community. 

The CKM-matrix element \V u b\ can be ex- 
tracted by comparing a measurement of the par- 
tial rate T U (P + < Ap) with a theoretical predic- 
tion for the product of the event fraction Fp(Ap) 
and the total inclusive B —> X u l~v rate. The re- 
sulting theoretical uncertainty on \V u b\ is 

^^ = (±7±4±5±4)%, (9) 



where the last error comes from the uncertainty 
in the total rate [13,14]. Because of the large 
efficiency of the P + cut, weak annihilation effects 
[11] have an influence on \ V u b \ of less than 2% and 
can be safely neglected. 
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